Spin-state transition and spin-polaron physics in cobalt oxide perovskites: 
ab initio approach based on quantum chemical methods 
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A fully ab initio scheme based on quantum chemical wavefunction methods is used to investigate 
the correlated multiorbital electronic structure of a 3d-metal compound, LaCoOs. The strong 
short-range electron correlations, involving both Co and O orbitals, are treated by multireference 
techniques. The use of effective parameters like the Hubbard U and interorbital [/', J terms and the 
problems associated with their explicit calculation are avoided with this approach. We provide new 
insight into the spin-state transition at about 90 K and the nature of charge carriers in the doped 
material. Our results indicate the formation of a t^gC^ high-spin state in LaCoOs for T > 90K. 
Additionally, we explain the paramagnetic phase in the low-temperature lightly doped compound 
through the formation of Zhang-Rice-like O hole states and ferromagnetic clusters. 



An accurate description of correlated electrons is one 
of the central problems of condensed-matter theory. The 
commonly used approach for the computation of elec- 
tronic structures is based on density-functional theory 
and various approximations to it, the most applied be- 
ing the local density approximation (LDA). This ap- 
proach has led to unprecedented progress in computa- 
tional condensed-matter physics and surprisingly good 
results for numerous compounds. However, it is also 
known and well documented that the LDA has prob- 
lems with describing strongly correlated systems which 
usually involve d and / electrons. For that reason ex- 
tensions like the LDA-|-?7 [l| and LDA plus dynamical 
mean-field theory (LDA-I-DMFT) ^ were developed. In 
the latter approaches, the strong local correlations are 
described by effective interactions such as the Hubbard 
U, an interorbital Coulomb repulsion [/', and an on-site 
exchange coupling J. Even so, the correlation treatment 
is limited to on-site effects, despite of the fact that inter- 
site correlations are known to be sizeable. To deal with 
the latter, the LDA+DMFT scheme has been generalized 
to a cluster LDA+DMFT 0, [3|. But a self-consistent 
computation of the Coulomb repulsion parameters U, U' 
(and intersite V) is in this context particularly difficult 
[j| . It is fair to state that the main virtue of LDA, namely 
its simplicity, has got more and more lost. 

This suggests the use of a different approach, based 
on the computation of many-body wavefunctions [5| by 
means of quantum chemical techniques [6| . It has the ad- 
vantage that all approximations are well controlled and 
that it is not necessary to introduce effective interactions 
like U, U' or J. Here, we demonstrate this for the case 
of a 3d oxide compound, LaCoOs. Our analysis is based 
on a local Hamiltonian approach [71, |8|, M UM ^^^d per- 
formed on large fragments including several CoOg oc- 
tahcdra. Localized Wannier functions and embedding 
potentials are obtained for these fragments from prior 
periodic Hartree-Fock (HF) calculations for the infinite 



crystal lOlilUl- Since our scheme is parameter-free, vari- 
ational, and multiconfigurational (or multireference, see 
below), it has great predictive power. It can, for example, 
predict the basic ingredients required in effective impu- 
rity models like DMFT. For late transition-metal (TM) 
oxides, in particular, it is important to anticipate whether 
(and which of) the TM d orbitals are sufficient for con- 
structing a minimal effective orbital space or whether 
both TM d and oxygen p functions must be considered 
as active. That such questions can be reliably addressed 
with quantum chemical methods was shown before for 
the layered copper oxides 1^ and ladder vanadates [l3| . 
However, simple point-charge embeddings were used in 
previous work for representing the surroundings of the 
region where the correlation treatment is carried out. A 
newly developed, rigorous embedding technique id. Ill| 
is applied here for the first time to a Sd-metal compound. 

LaCoOa has attracted considerable attention due to a 
number of puzzling phase transitions induced by changes 
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m temperature 

18|. Up to now 



15j . doping 15J, Il6|, and/or strain 
most of the experimental work 



was aimed at understanding the nature of the phase 
transitions in the undoped compound, from a nonmag- 
netic insulator at low T to a paramagnetic semiconduc- 
tor above 90 K and to a metal for T> 500 K. The low-T 
ground-state was assigned to a closed-shell ^20^" config- 
uration [low-spin (LS), 5 = 0] of the Co ions [14,]. For 
T > 90 K, however, the available experimental results 
are rather contradictory. While recent x-ray absorption 
spectroscopy (XAS) [19| and inelastic neutron scatter- 
ing (INS) [20| measurements indicate with increasing T 
a gradual transition into a 5 = 2 (f^ge^) high-spin (HS) 
configuration of the Co 3d electrons, electron energy-loss 
spectroscopy data 1211 and the observation of Co-0 bond- 
length alternation [22| suggest the formation of a S" = 1 
(tlgCg) intermediate-spin (IS) state for T>90 K, suscep- 
tible to Jahn- Teller distortions. On the theoretical side, 
the transition to a HS state around 90 K is favored by 



many-body model Hamiltonian calculations [l9|. The IS 
electron configuration is supported by LT)A+U calcula- 
tions [2l|, [23] ■ A wavefunction-based quantum chemical 
analysis is not available yet for this system. 

The first step in our investigation is an ab initio closed- 
shell restricted HF (RHF) calculation for the periodic 
solid. The CRYSTAL code [ij] is employed for this pur- 
pose. We use all-electron Gaussian-type basis sets of 
triple-zeta quality 25| and a cubic perovskite crystal 
structure. The Co-0 distance is taken to be d = 1.92 
A, as reported for low T in ref. [28|. It is well-known 
that the HF approximation strongly overestimates the 
fundamental band gap in insulators. For LaCoOa, the 
RHF gap is 13.3 eV, where the upper valence bands have 
O 2p character and the lowest conduction bands are re- 
lated to the Co 3d-eg levels. On-site and intersite corre- 
lation effects are here investigated in direct space by mul- 
ticonfiguration complete-active-space self-consistent-field 
(CASSCF) and multireference configuration-interaction 
(MRCI) calculations p]. The orbital basis entering 
the correlation treatment is provided by the localization 
module of the crystal program. These functions are the 
Wannier orbitals (WO's) associated with the RHF bands 
[13 . The CASSCF and MRCI calculations are carried 
out with the molpro quantum chemical package [3C| . 

Several sites are included in the correlation treatment. 
These sites are partitioned into two groups: an "active" 
region Ca, which here consists of one CoOg octahedron, 
and a "buffer" zone Cb, including the Co and all other 
O ions of the six nearest-neighbor (NN) octahedra (i.e., 
6 Co and 30 O ions) plus the La NN's of the "active" 
octahedron (8 La sites) 3l| . The role of the buffer zone 
Cb is to ensure an accurate representation of the tails 
of the WO's centered in the active region Ca- While 
the occupied orbitals in the buffer zone are frozen, or- 
bitals centered at sites in the Ca region (and their tails 
in Cb) are allowed to relax and polarize in the quantum 
chemical calculation. We denote Ca+Cb as C. The sur- 
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FIG. 1; Diagram for the lowest A^ (left) and (iV-1) (right) 
particle states. The MRCI treatment is based on a cas-7 
reference, see text. Only the Co 3d and O "ip-Cg electrons are 
correlated in the MRCI-7 calculations. For MRCI-11, the Co 
3s, 3p electrons are also included. The reference is in each 
graph the energy of the MRCI-11 LS state. 



roundings of this group of ions (i.e., the crystalline envi- 
ronment) are modelled as an effective one-electron poten- 
tial. This is obtained from the periodic RHF calculation 
and incorporated in the one-electron Hamiltonian in the 
CASSCF/MRCI treatment via the CRYSTAL-molpro in- 



terface program [lO|, |ll| . A ground-state closed-shell self- 



consistent-field calculation for C leads to changes of only 
few meV (« 3 meV) in the total energy, which shows 
that our embedding potential was properly constructed 
and the projection of the original crystal WO's centered 
at sites of the Ca region onto the set of basis functions 
associated with C — Ca+ Cb |32| | leads to insignificant 
changes of the longer-range "tails" . 

The low-temperature (T ^ 90 K) ground-state config- 
uration and the nature of the lowest A^-particle excited 
states are first studied at the CASSCF level. In CASSCF, 
a number of active electrons is allowed to be distributed 
in all possible ways over a given number of active orbitals. 
Both the orbitals and the coefficients in the multiconfig- 
uration expansion are optimized [6|. In the case of a 
single Co^"*" ion, a minimal active space is constructed 
with six electrons and five 3d orbitals, which we denote 
as cas-5. Intriguingly, for such a minimal active space, 
CASSCF predicts a HS tjge'i {^T2g) ground-state config- 
uration. The LS state, t2g^g (^^3), is 1.26 eV higher in 
energy. The IS state, tlgc] {^Tig), is 0.21 eV above the 
LS state, see Fig. 1. Clearly, the minimal orbital space 
description is insufficient for this system. 

A first guess is that the active space must be enlarged 
with O 2p orbitals. Indeed, the RHF band structure, 
with the O 2p bands located above the Co t2g levels, 
see above, suggests that O 2p to Co 3d charge-transfer 
(CT) effects are important. The two a-like eg combina- 
tions of 2p orbitals at the NN O sites have the largest 
overlap with the Co 3d functions. When these two com- 
binations of Cg symmetry are added to the active space, 
which we denote as the cas-7 orbital space, the CASSCF 
HS-LS splitting is reduced by 0.39 eV, to 0.87 eV. This is 
related to the strong interaction between the "non-CT" 
tjgCg and the O 2p-eg to Co 3d-eg CT configurations in 
the LS CASSCF wavefunction. However, the HS state is 
still the lowest, see Fig. 1. Other types of charge fluctu- 
ations are taken into account by MRCI calculations. In 
a first step, we account for all single and double excita- 
tions from the Co 3d and O 2p-eg orbitals. The reference 
is the cas-7 wavefunction. The HS-LS splitting is now 
reduced by an additional 0.69 eV, with the HS state re- 
maining the ground-state. A switch in the energy order 
of these two states is only obtained when accounting for 
correlation effects (single and double excitations in the 
MRCI treatment) which involve the semi-core Co 3s, 3p 
electrons. These correlations are mainly related to the 
coupling between the O 2p-eg to Co Sd-Cg CT, on-site 
Co 3s, 3p -^ 3d, 4s,4p, and Co 3d -^ 4s, 4p excitations. 
For such MRCI wavefunctions [34] , the LS state is found 
indeed to be the lowest and the first iV-particle excited 



state is the HS state, with an excitation energy of 6 meV 
(«70 K, not visible in Fig. 1). 

We thus find that the HS state is lower in energy than 
the IS arrangement. For an undistorted lattice, this 
energy diS'erence is 0.60 eV. Additionaly, the quantum 
chemical calculations provide unique, detailed informa- 
tion about the importance of different types of correla- 
tion effects. The good agreement between the LS-HS 
splitting and the temperature where the magnetic phase 
transition occurs (70 vs. 90-100 K) is remarkable. One 
should, however, realize that inclusion of other excitation 
processes and larger basis sets would modify quantita- 
tively this splitting, although we expect that the order of 
the two states would not change [35J . 

We now proceed with the analysis of the lowest 
electron-removal and electron-addition states. It was ar- 
gued that the largest interaction with the 2p -^ 3d CT 
configurations and the strongest stabilization effects oc- 
cur for the ^20^0 IS Iiols state [3a]. Therefore, the lowest 
(A'^— 1) states are often associated with an IS, f^gSg {"^Tig), 
configuration of the Co ions. Our calculations show that 
the lowest-energy hole states are related to the LS, ^L^o 
{^T2g), configuration. Even by CASSCF, with a cas-7 
active space, the LS ^io^o ^ol^ state is about half eV 
lower than the IS {N—1) state. MRCI calculations which 
include the Co 3s, 3p electrons yield an even a larger 
splitting of 0.85 eV, see Fig. 1. The implication is that 
the strongest CT effects occur for reference configura- 
tions with completely empty 3d-eg levels, as also found 
in the iV-electron system. 

The lowest hole states in LaCoOa have strong O 
2p character and resemble the lowest {N — 1) state 



in cuprates 37|. In cuprates, the dominant contri- 



bution to the ionized state is a superposition of two 
configurations, \t^gdl2b'^a'^) and jf^gd^afe^a^), where b 
and a are bonding and antibonding combinations of 
the O 2po- and Cu 3d orbitals of {x^ — y^) symme- 
try on one CUO4 plaquette. In the CASSCF wave- 
function 12], the Cu "id and O 2pcr (x'^—y^) orbitals 
contribute with nearly equal weight to b and a. The 
added hole has therefore predominant O 2p character 
[13]. For cubic perovskites, the two id-Cg components 
are degenerate. The dominant contributions to the first 
ionized state in LaCoOa imply then three configura- 

z^ x^—y^z^ I ^ I ^g x^—y^z^ x-^—y^z^l^ 



tions: \t\gb\i. 
and lt2g^^2_j,2 



'^2_„2a?2)- The bonding and antibond- 

X y z I 

ing Cg orbital combinations are again strong mixtures of 
TM 3d and O 2p functions (see Fig. 2), such that the 
added hole has large weight at the O sites. However, 
in contrast to the hole state in cuprates, all six ligands 
are involved, although in the presence of lattice distor- 
tions [22] the eg levels split and the hole does not have 
the same weight at each of the NN ligands. A study of 
the coupling between the charge and lattice degrees of 
freedom in LaCoOa is left, however, for future work. 
The strong O 2p hole character of the ionized state is 



expected to induce ferromagnetic (FM) correlations be- 
tween the Co i^-like ion and its six Co NN's, like for the 
2p hole in cuprates [38| or for the O 2pJ configuration in 
NaV205 [l3|. The formation of FM spin-polaron clusters 
was indeed inferred from measurements of the magneti- 
zation in lightly doped samples, see, e.g., ref. :15:]. For 
no or small distortions, the configuration of each of the 
six Co NN's is HS ^\ge^g and the total spin of the spin 
cluster is 5 = 1/2 + 6 x 2 = 25/2. If the coupling to the 
lattice is strong, an IS configuration t\ge^g at the NN Co 
sites (22| can be energetically more favorable and the to- 
tal spin is 5=1/2+6x1 = 13/2. Obviously, the motion of 
the added hole is strongly renormalized by the NN FM 
correlations. In cuprates, for example, the NN "bare" 
hopping is reduced by a factor of 4 by short-range spin 
interactions I4] . The study of such correlations require 



very involved calculations for LaCoOs because a number 
of five d orbitals is needed at each TM site. A computa- 
tional approach like that used in copper oxides, where all 
TM neighbors of the two octahedra directly involved in 
the hopping process are included in CASSCF [l2], is at 
present not feasible for LaCoOs. The possibility of using 
an incremental scheme like that proposed in ref. [7,] will 
be investigated in future work. 

The lowest (A^+1) states are related to the HS ^20^0 
(■^Tig) configuration. MRCI calculations which include 
in addition to the orbitals used in the cas-7 reference the 
Co 3s, 3p predict a splitting of 1.00 eV between the t\ge^g 
{^Tig) and tl^ej {^Eg) states. 

The matrix elements (ME's) associated with local 
electron-removal and electron-addition configurations 
like those described above can be used to construct "cor- 
related" band structures within a quasiparticle picture 
[34]. The N -^ (A^+1) excitation energies correspond to 
the diagonal ("on-site") ME's in the k-dependent secu- 
lar problem while overlap and Hamiltonian integrals in- 





b) 





FIG. 2: Bonding and antibonding p-d orbitals for the LS hole 
state of (a) (a;^ — y^) and (b) z^ symmetry. The (x'^ — y^) 
hybrids [see (a)] resemble the ZR p-d orbitals in cuprates. 



volving holes or added electrons at different lattice posi- 
tions R (i.e., (^l^^^l^^J^J and (*^„=f;|^|^^t1j mE's, 
where n is a band index and a denotes the spin) enter 
as ofF-diagonal (intersite) terms. As discussed above, ac- 
curate results for the intersite ME's cannot be obtained 
at this moment in LaCoOs. Nevertheless, we can extract 
from our data estimates for the ionization potential [low- 
est A^ — > (TV— 1) excitation energy] and electron affinity 
[lowest (AM-1) state]. In our MRCI calculations, these val- 
ues are IPq — 4.47 and EAq — —4.24 eV and incorporate 
the most important short-range correlations. Corrections 
due to long-range polarization effects can be estimated 
by applying the approximation of a dielectric continuum 
[ai]. This gives IP = IPo - CIR and EA = EAq + C/R, 
where C = e^^^^, eo ~ 23 [3§] is the static dielectric 
constant, and R defines a sphere around the extra hole 
or electron beyond which the dielectric response reaches 
its asymptotic value eg- R is taken here as the average 
between the NN Co-0 and Co-La distances, i? = 2.62 A, 
which leads to IP = 1.84 and EA = -1.61. Given the 
fact that we neglect the finite widths of the two bands, 
the energy difference IP — EA = 3.45 eV represents an 
overestimate of the actual band gap. 

To summarize, a recently developed ah initio scheme 
is applied here to the study of many-body effects in 
LaCoOa. The nonmagnetic low-T ground-state involves 
a subtle interplay between ligand-field and exchange in- 
teractions plus on-site and CT fluctuations. The lowest 
A^-particle excited state is predicted to be the HS t'^gei 
state. In agreement with XAS [19] and INS [gfl] exper- 
iments, the IS i2gCg state is 0.60 eV higher. Whether 
this HS-IS splitting is sufficiently low to be compensated 
by the energy gain associated with a Jahn- Teller ordered 
configuration of the Cg orbitals, as suggested in ref. [22], 
will be the subject of future work. We also study the 
nature of the low-energy electron-removal and electron- 
addition states. The lowest (A^ — 1) excitation is related 
to a local tig J '5'= 1/2 configuration. The added hole has 
large weight at the O sites, due to a strong mixing of 
the Co and O eg orbitals. A large hole density at the O 
sites causes FM correlations among the adjacent Co Zd 
electrons, which agrees with the observation of "giant" 
FM polarons (5 = 10—16) a,nd paramagnetic behavior at 
small doping and low T [15| . Our analysis offers new in- 
sight into the correlated electronic structure of LaCoOs. 
Additionally, the ah initio results provide valuable infor- 
mation for models aiming at a realistic description of the 
finite-temperature properties of this system. 
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